Anisotropic superconductivity in the two-dimensional electron-phonon system 
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The properties of the d-wave superconducting state in the two-dimensional system have been 
studied. It has been assumed, that the pairing mechanism is based on the electron-phonon and the 
electron-electron-phonon interactions. The obtained results have shown the energy gap amplitude 
(Atot) crossover, from the BCS to non-BCS behavior, as the value of the electron-electron-phonon 
potential increases. The model has been tested for the La2-xSrxCu04 and Bi2Sr2CaCu208+i high- 
Tc superconductors. It has been shown, that the dependence of the 2Aj^j/fc_BTc ratio on the hole 
density is in agreement with the experimental data. 

PACS numbers: 74.20.-z, 74.20.Fg, 74.20.Mn, 74.20.Rp, 74.25. Bt, 74.72.-h 



In year 1986, Bednorz and Miiller have discovered the 
high temperature superconductivity [ij. Despite of the 
large efforts, the origin of the pairing correlations in 
high-Tc materials is still unclear. In particular, the de- 
scription of the cuprates' properties in the framework 
of the pure electronic models (the Hubbard or related 
approaches) is probably insufficient, since in these mod- 
els the pairing correlations are too small On the 
other hand, the conventional electron-phonon interaction 
is weak enough, that the depairing electronic correlations 
should completely suppress the phonon-mediated super- 
conductivity [3]. 

In order to solve the discussed problem, we assume, 
that the pairing mechanism would inseparably link to- 
gether the strong electron correlations and the crystal lat- 
tice vibrations (the electron-phonon (EPH) and electron- 
electron-phonon (EEPH) interaction). We notice, that 
the above idea has been originally postulated in Q , where 
the properties of the anomalous s-wave superconduct- 
ing state have been studied. The obtained results have 
shown, that for the high value of the EEPH coupling, 
the s-wave energy gap weakly depends on the temper- 
ature (T < Tc), and above the critical temperature it 
vanishes at the Nernst temperature. 

In the presented paper we will study the thermody- 
namic properties of the d-wave superconducting state on 
the basis of the postulated pairing model. In particu- 
lar, we will calculate the dependence of the energy gap 
amplitude Atot on the temperature for the selected val- 
ues of the EEPH potential. Next, we will analyze the 
dependence of the 2A[^l/kBTc ratio on the hole den- 
sity for La2-xSrxCu04 (LSCO) and Bi2Sr2CaCu208+5 
(Bi2212) superconductors. The theoretical predictions 
will be compared with the experimental data. 

The effective Hamiltonian takes the form Q: 



(1) 



The first term represents the non-interacting electrons: 
i7(°)^5]£kcLck., (2) 

kcr 

where: Sk = —tj (k) denotes the electron band 
energy for the two-dimensional square lattice; t is 
the nearest-neighbor hopping integral and 7 (k) = 
2 [cos (kx) + cos {ky}]. The symbol c|^^, (cko-) is the cre- 
ation (annihilation) operator for the electron with mo- 
mentum k and spin a. The EPH and EEPH interaction 
terms are given by: 
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and 
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functions Vkq and t/kk'qi indicate the pairing potentials: 
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where g'-^^ and g'^^^ represents the EPH and EEPH cou- 
pling respectively; ujq denotes the characteristic phonon 
frequency. 

On the basis of the operators ^ and ^ it is pos- 
sible to deduce the Hamiltonians which describe the d- 
wave superconducting state. In the case of the Hamil- 
tonian ([3]) this procedure is known and widely described 
in the literature (see e.g. the paper 0). With reference 
to the above, we will discuss only the derivation of the 
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Now, we assume: U, 



(-) 



We 



notice that, the number 24 has been introduced, since 
the potential energy term represents the interaction be- 
tween every four of particles counted once; Nq is the 
normalization factor: Nq = l/X^k"' '^here the symbol 
Yl^' denotes the sum over the states for which |ek| < i^o 
(the BCS-type cut-off). In the next step, we limit the 
symmetry of the energy gap to the dominating d-wave 
symmetry. The total Hamiltonian after applying the ap- 
proximation presented in the paper Q takes the form: 



kcr 



(11) 



E[^k''^4tc'-k, + Af)c_k,Ckt 



FIG. 1: The dependence of the energy gap amplitude on the 
temperature for the selected values of the EEPH potential. 
The solid line represents the physical stable solution; the dot- 
ted line corresponds to the unstable solution, where the ther- 
modynamic potential is bigger than in the first case. The 
vertical line indicates the position of the critical temperature. 



anisotropic Hamiltonian on the basis of the EEPH oper- 
ator. In the first step, we separate the momentums in 
the expression 



The symbol F^'') and U^'^'> denotes the d-wave effec- 
tive potential for the EPH and EEPH channel re- 
spectively. In particular: V'^'^^ = V/2 and U'^^^ = 
U/8. The anisotropic order parameter is given by: 



A 



= a'''-*?] (k), where the amplitude is expressed 
as: A('') s i^Er^(k)(c-k;Ckt) and (k) ee 
2 [cos (kx) — cos {ky)]. 

On the basis of the Hamiltonian ([TT|) , we calculate the 
thermodynamic Green function by using the equation of 
motion method ra . The result has the form: 



ki ~k4(T 
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(7) 

we transform the where 



where: /ikgkicr = ^a-o-^-kso-'^-kicrCki-cr. With help of 
the relation: Ck<T = J2j e^^^'^'Cj^, 
operator ([7]) to the Wannier representation, where we re- 
strict ourselves to on-site and the nearest neighbor pair- 
ing. The Hamiltonian takes the form: 
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where h 
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and h 



(+) 

( + ) 



c;3_^Cj4-ac;-^^Cj-3^. The symbols U^^^^ and U^^^^ de- 
note the local and kinetic potential respectively. Next, 
we return to the Bloch representation: 
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We turn the reader's attention toward the fact, that the 
obtained Green function possesses the analytical struc- 
ture, which is more complex than the structure of the 
BCS Green function [1], Q- In particular, the energy 
gap amplitude is the complicated function of the order 
parameter amplitude. However, the energy gap, in spite 
of its complicated form, is characterized by the pure d- 
wave symmetry. 

On the basis of Eq. p2|) we derive the fundamental 
thermodynamic equation: 
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The sum X]k° approximated in the following manner: 
Lk" - /-TT /-TT dkxdkyO {ujQ - Isk^^kyl), where the sym- 
bol 9 represents the unit step function. In the model 
calculations, we have taken t as the energy unit. 

In Fig.[T] we present the temperature depen- 
dence of the energy gap amplitude (Afot = 
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TABLE I: The parameters V'^'^^ and (7^''' calculated by using Tc and the mean values of T*. 



Material 


Type 


t (meV) 


Ref. 


ojo (meV) 


Ref. 


Tc (K) 


T* (K) 


Ref. 


l/'"^ (meV) 


U"*"' (meV) 
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p=0.15 










38 


135.5 ±30.2 


[10] 


4.41 


53.89 




p=0.22 










28 


44.3 ± 30.2 


[10] 


3.91 


45.11 


Bi2212 


p = 0.125" 


350 


[12], [13], 


80 


[15], [16], 


83 


~ 290 


[20] 


5.47 


55.40 




p = 0.143 




[14J. 






90 


~ 300 


[21] 


5.68 


56.37 




p = 0.193 








fl9]. 


84 


~ 270 


[21] 


5.51 


51.72 




p = 0.198 










81 


~ 220 


[21] 


5.40 


43.88 



"The hole density p has been estimated as the doping p = x. 

''The hole density p has been obtained on the basis of the empirical 
formula Tc (p) /Tc,max = 1 - 82.6 (p - 0.16)^ where Tcmax = 
92.2 K 
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and the selected values of JJ*^''^. It is easy to see, that 
for the high values of the EEPH potential, the shape of 
the function Atot (T) is sharply different from the BCS 
prediction. In particular, for T S (0, Tc) the energy gap 
is very weakly dependent on the temperature; up to the 
critical temperature Atot extends into the anomalous 
normal state to the temperature T* . In the case of 
the c?-wave superconducting state, the temperature 
T* is interpreted as the pseudogap temperature (in 
contrast to the s-wave superconductivity, where the 
highest value of the temperature, for which the non-zero 
solution of the gap equation exists, is connected with the 
Nernst temperature T** We notice, that usually: 

T^ ^ 

Below we compare the theoretical predictions with the 
experimental data for LSCO and Bi2212 superconduc- 
tors. For this purpose, we have calculated the values of 
the pairing potentials on the basis of Tc and T* experi- 
mental values. The obtained results have been collected 
in Tab.lll Next, by using the T^^'') and U'-'^'^ values, the 
hole density dependence of the Ri = 2A[^l/kBTc ratio 
has been obtained. We notice, that the energy gap am- 
plitude at the temperature of zero Kelvin is defined as: 



,(0) 



6 



Av) 



where Aq'''' denotes 



the order parameter amplitude at the temperature of zero 
Kelvin. 

In Fig. [5] we present the dependence of the Ri ratio 
on the hole density for LSCO superconductor. It can be 
seen, that with the increase of p, the parameter Ri suc- 
cessively decreases. In particular, for the underdoped re- 
gion {p < 0.155) the values of Ri are significantly higher 
than the d-wave BCS value 4.28 [H]. Slightly above 
p = 0.155 (the overdoped region) the Ri ratio approaches 
closely the weak-coupling c?-wave BCS result. We notice, 
that for p > 0.15, the lower accuracy of the theoreti- 
cal results can not be determined, since Tc and T* are 
experimentally indistinguishable (see Tab. IJ). Now, we 
have compared the theoretical predictions which the ex- 
perimental values of -Ri, received by the few different 
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FIG. 2: The ratio Ri as a function of p for LCSO supercon- 
ductor. The solid lines with the open squares represent the 
theoretical calculations based on the data presented in the 
paper [l^. The overshadowed areas mean the accuracy of 
the achieved results. The filled symbols correspond to the ex- 
perimental results obtained by: (a) - Hashimoto, et al. p^ . 
(b) - Nakano, et al. (c) - Oda, et al. (d) - Kato, et 

al. [23], (e) - Wang, et al. fSslj, (f) - Yoshida, et al. [H], (g) 
- Wen, et al. 



researchers. The qualitative agreement of the theoretical 
predictions with the experimental data proves, that the 
measured dependence of the ratio i?i on p can be well 
reproduced with an use of the presented model. 

In Fig. [3] we show the shape of the function Ri (p) 
for Bi2212 superconductor. The presented results prove, 
that the theoretical line determines the high value of 
Ri in the whole range of the considered hole density; 
p e (0.125, 198). Important is also the fact, that the 
model correctly reconstructs the experimental data. 

To summarize, we have presented the model that de- 
scribes the properties of the d-wave superconducting 
state in the two-dimensional system. We have assumed, 
that the pairing correlations are induced by the electron- 
phonon and electron-electron-phonon interactions. In the 
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FIG. 3: The dependence of the ratio Ri on p for Bi2212. 
The sohd hne with the open squares represents the theoretical 
calculation. The filled and half-filled symbols correspond to 
the experimental results obtained by: (a) - Renner, et al. [20j . 
(b) - Hoffmann, et al. .29] , (c) - Ponomarev, et al. [13], (d) - 
Oki, et al. [sj], (e) - Krasnov, et al. [s^], (f) - Gupta, et al. 
[3^ . (g) - Kanieel, et al. [sj, (h) - Carnpuzano, et al. [ssj ]. 
Tanaka, et al. [3a ]. (i) - Nakano, et al. 22|, (j) - Oda, et al. 
[13], (k) - McElroy^i al. [sl], (1) - Matsuda, et al. [sg, (m) 
- Hoffman, et al. [40|. The lines (n) were obtained by using 
the empirical relation: Ri (p) = (15 ± 1) - (38 ±5)p [4l|]. 



first step, we have derived the fundamental thermody- 
namic equation. Next, on the basis of the exact numeri- 
cal solution, we have shown, that for the high value of the 
EEPH potential, the temperature dependence of the en- 
ergy gap amplitude differs sharply from the BCS predic- 
tion. In particular, the energy gap amplitude is slightly 
dependent on the temperature for T g (0, Tc); above the 
critical temperature, the energy gap amplitude persists to 
the pseudogap temperature. In the paper the theoretical 
predictions have been compared with the experimental 
data for LCSO and Bi2212 superconductors. It has been 
shown, that the calculated hole density dependence of the 
i?i ratio correctly reproduces the experimental results. 
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